Induced P-wave Superfluidity in Asymmetric Fermi Gases 
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We show that two new intra-species P-wave superfluid phases appear in two-component asymmet- 
ric Fermi systems with short-range S-wave interactions. In the BEC limit, phonons of the molecular 
BEC induce P-wave superfluidity in the excess fermions. In the BCS limit, density fluctuations 
induce P-wave superfluidity in both the majority and the minority species. These phases may be 
realized in experiments with spin-polarized Fermi gases. 
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The phenomenon of fermionic superfluidity spans 
twenty orders of magnitude, from cold atomic gases, 
through liquid ^He, electronic superconductors, nuclei, 

' and neutron stars, to quark color superconductors. Ex- 
periments with trapped cold atoms provide clean access 

\ to this physics, with tunable interactions and composi- 
tions. 

In this Letter, we discuss the possibility of P-wave 
superfluidity in cold asymmetric Fermi gases compris- 
ing two species (t and i) of equal mass m, with short- 
range S-wave interactions. (The generalization to un- 

\ equal masses is straightforward.) If the system is sym- 
metric (equal number density of each species) , the ground 

\ state is well described as a fully gapped superfluid: There 
is a crossover from a BCS state at weak attraction to 
a BEC of tightly-bound diatomic molecules (dimers) at 
strong attraction. 

\ The ground state of asymmetric systems is not even 
qualitatively well understood. If the chemical potential 

' difference 5/i between the two species is larger than the 
gap A in the spectrum, the stability of the superfluid 

\ state is compromised. At weak coupling it is known 
that there must be a phase transition for (5/i ^ A/\/2 

] to some asymmetric state Various possible asym- 
metric superfluid phases have been proposed: For ex- 
ample, anisotropic/inhomogeneous superfluid states with 
crystalline structure (LOFF) 5| or deformed Fermi sur- 
faces yL and homogeneous gapless (breached pair) super- 
fluids y|. To date, all proposed phase diagrams conclude 

\ that, for large asymmetry, one finds a normal Fermi liq- 

■ uid Bull 

\ Here we show that attractive mtra-species P-wave in- 
teractions are induced between gapless fermions, leading 
to the formation of P-wave superfluids. We propose a 
phase diagram (Fig.0 where phases with excess fermions 
are P-wave superfluids (except for the free Fermi gas con- 
sisting of a single non- interacting species). In the BEC 
regime, the Pi phase consists of a molecular BEC with 
a single P-wave superfluid of the excess particles. In the 
BCS regime, the P2 phase has two coexisting P-wave su- 
perfluids: one in the majority component, and one in the 
minority component. 

We calculate the induced interactions that give rise to 
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FIG. 1; Conjectured zero-temperature grand canonical phase 
diagram after 6] . Scaling relations allow this two-dimensional 
projection of the three parameters {jj,,a,5fj,). The chemi- 
cal potential asymmetry Sfi — {/i^ — /i^)/2, and the S-wave 
scattering-length a are expressed in units of the S-wave gap 
A^fi^o and density na^^o of the symmetric state with the 
same parameters /i and a, but with Sfi = 0. The phase 
Pi — a gapless superfluid with single Fermi surface in |g — is 
our P-wave superfluid/BEC phase. The phase P2 — a two- 
component Fermi-liquid in (Qj — is our two-component P-wave 
phase. There may be additional P-wave phases (for instance 
in the breached pair phases) near the S-wave Feshbach res- 
onance. Likewise, the exact nature of the phase(s) labeled 
LOFF? has yet to be determined. This phase diagram is 
illustrative: Our results are independent of its quantitative 
structure. The dotted lines in the BEC and BCS regimes 
represent sample trajectories of constant a and Sjj,. A trap 
provides a radially varying fi, and would comprise concentric 
shells of the various phases crossed by these trajectories start- 
ing from the crosses at the center of the trap. To maximize 
the physical space occupied by these new phases, it might be 
useful to load the trap so the core is a P-wave phase as indi- 
cated by the two shorter trajectories. (Note: we have verified 
the qualitative structure and shape of this diagram (outside 
of the unitary regime) with a self-consistent mean-field cross- 
over model.) 



P-wave superfluidity in the controlled limit of small scat- 
tering length, and estimate the dependence of the P-wave 
gap Ap on the various physical parameters. These P- 
wave phases could soon be observed, as experiments are 
just starting to explore the properties of spin-polarized 
gases II, 0- 
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BEC Regime. In the BEC limit of strong attrac- 
tion, the fermions form tightly-bound molecules that con- 
dense into a BEC. For asymmetric gases, the naive 
picture is that the excess fermions interact weakly with 
the molecules, and coexist in the "spaces" between the 
molecules. The resulting phase exhibits gapless fermionic 
excitations in a superfluid BEC and has been called 
a "gapless superfluid". (This picture is supported by 
Ref. mj and mean-field cross-over models.) We show, 
however, that phonon exchange induces an attractive P- 
wave interaction between the excess fermions. This in- 
teraction destabilizes the Fermi surface, forming an in- 
fra-species P-wave condensate on top of the molecular 
BEC. 

To estimate the magnitude of the induced P-wave in- 
teraction Up, and the relevant critical temperature T^f , 
we consider the exchange of a sing le p honon between the 
excess fermions. In weak coupling |2a|, the P-wave super- 
fluid pairing gap is given by Ap ~ ep / (NpUp)] , 
where ep is the Fermi energy of the excess fermions, 
A'p — mkp / {2Tr'^h?) is the density of states at the Fermi 
surface, and hkp is the Fermi momentum. The interac- 
tion induced by single phonon exchange is given by \vi\ 



UindiqoA) = Ul 



(1) 



where t/fb is the fermion-boson coupling, Uhh is the 
boson-boson coupling, E^ = g^/(2mb), ?7ib = 2rn is 
the mass of the molecules, and rib is the boson density. 
The couplings are determined from the scattering length 
a > using few-body methods [l^: 
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We neglect the frequency dependence in weak coupling, 
and thus the static induced interaction is 
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where the speed of sound in the Bose gas is 



c = y/ubUhh/ rrib = hy/2T:a Tib abb / ^b ■ 



(3) 



(4) 



In weak coupling, the infrared behavior implies that mo- 
menta in the vicinity of the Fermi surface |pi| = hkp 
dominate the pairing interaction. Moreover, the en- 
ergy transfer is small, justifying our neglect of the fre- 
quency dependence. Pauli exclusion forbids S-wave in- 
duced interactions between the excess fermions, thus, 
the induced P-wave interaction is dominant. Project- 
ing onto the Legendre polynomial Pi (cos 0) with cos 9 — 



1 — (pi — p2)^/(2fc|,) leads to the P-wave component 



Up = 



dcos( 



cos 6* C/ind(0,Pl - P2) : 



m abb 



mbC 



where 



2 








{x^ " 


^0 



Ri{x) 

The asymmetry enters through the parameter 
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where np — fcp/(67r^) is the density of excess fermions. 
The function Ri{x) has a maximum at i?(1.86) = 0.104 
and the following limiting behavior: 



Ri{x) 
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ln(2;^)/a;^ where a; 3> 1 



where a; <C 1 , 
2 1 



For small asymmetries, np/n^, <^ kpa <C 1, we find the 
intriguing result that the P-wave gap is independent of 
the scattering length, 
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although it is suppressed for small asymmetry. For large 
asymmetries, np/ny, ^ kpa, the P-wave gap is 
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If we fix the interaction and density of excess fermions 
(constant a and kp), we find a maximal P-wave gap as a 
function of the boson density: 
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for a small intermediate asymmetry: 

Tip /lib ~ 0.44 kpa < 1 . 
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Note that this maximal P-wave gap depends parametri- 
cally on l/(fcpa) in the exponent as for S-wave pairing 
in the BCS regime. Therefore, we expect similar high 
critical temperatures towards resonance. 

It is instructive to compare the critical temperature 
for P-wave superfiuidity ~ Ap to the critical tem- 
perature for Bose condensation of the molecules T^^^ ~ 
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n^''^/(C(3/2)^/^TOb)- For the maximum P-wave gap, we 



find, with Eq. l(TT|l and (fT^ . the ratio 

rpP, 



J^BEC 
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which increases toward resonance. (Neglecting higher or- 
der effects, the prefactor may be calculated and is of or- 
der unity.) We therefore conclude that the P-wave su- 
perfluid/BEC phase Pi may be observed in asymmetric 
Fermi gases as one approaches the S-wave Feshbach res- 
onance from the BEC regime. 

Our results are justified in weak coupling. The omitted 
pre-exponential factors, however, depend on higher-order 
induced interactions. It is expected that the resulting 
factors are of order unity. We also note that the pairing 
interaction induced by single phonon exchange is related 
to the effective mass m* /m = 1 -I- Fi/3 through the Lan- 
dau parameter Fi = —NfUp. Effective mass corrections 
are thus higher order, but they would increase the den- 
sity of states at the Fermi surface, and thus increase the 
magnitude of the P-wave gap and . 

BCS Regime. Next, we show that all proposed asym- 
metric Fermi liquid phases are unstable towards a two 
component P-wave superfluid due to the exchange of den- 
sity fluctuations. This occurs in the BCS regime and is 
denoted by P2 in our phase diagram Fig. We start 
from a two component asymmetric Fermi gas with Fermi 

T I 

momenta fcp > fcp, and calculate the induced interactions 
in weak coupling. To lowest order in the S-wave interac- 
tion ATrah?/m, the induced interaction for back-to-back 
scattering is given by |2^ : 
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P-wave pairing in neutron stars, it is known that cen- 
tral induced interactions are attractive , but repulsive 
spin-orbit fluctuations dominate this effect [l^ . 

The resulting P-wave superfluid gap for the majority 
component is given by Ap ~ exp[l / {NpUp )] , where 
up denotes the P-wave projection of the induced inter- 
action as in Eq. ^ '2T|. This leads to 
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with the P-wave superfluid gap for the minority compo- 
nent given by interchanging the spin labels. The asym- 
metry enters through the function 



which has the limiting behavior 
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{z'^/18 where z <C 1 , 

K-F (7-41n2)(z- 1)/15 where z « 1 , 
21n(z)/(3z2) where z>l. 

For the symmetric case, we recover the result of 15], 
Ap £f exp[-7rV(4fc|a2K)], with k ^ (21n2 - l)/5, 
however, the work of 'l3| considered repulsive S-wave in- 
teractions. In our case, inter-species S-wave pairing will 
dominate for the symmetric system. The phase P2 will 
start for some small but finite asymmetry, and the devi- 
ations in the exponent will be linear in (z — 1). 

For large asymmetries fcp 3> fcp, the P-wave gap of the 
majority component is 



exp 
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while that of the minority component is 



(17) 



L(|pi-p2|/(2n4)). (14) 



The induced interaction for the minority fermions is ob- 
tained by interchanging the spin labels. As before, in 
weak coupling we neglect the frequency dependence and 
consider momenta on the Fermi surface. Thus, L{y) de- 
notes the static Lindhard function 
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The importance of induced interactions for superfluidity 
has been pointed out for symmetric Fermi systems: in 
weak coupling for S-wave pairing [l^ , for P-wave pairing 
with repulsive interaction and close to the Fesh- 

bach resonance (based on 0|). In addition, it has 
been shown that induced interactions significantly sup- 
press the superfluid gaps in neutron stars [l^ Il9| | . For 
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The majority component has a larger gap, but both are 

suppressed for large asymmetry. 

For fixed fcp, the minority gap Ap decreases monoton- 

t 

ically for increasing asymmetry, while for fixed fcp, the 
majority gap Ap has a maximum at fcp 



Vp has a maximum at fcp 
to the maximum of Li(z)/z = 0.11 for z ~ 1.3: 



0.77 fc][,, due 



exp 



0.11 (2fcjt,a) 



(19) 



Finally, we note that the P2 phase does not destabilize 
LOFF, or similar phases, whose condensation energy is 
parametrically the same as that of the S-wave BCS phase 
where As ~ exp(7r/2 fcpa). Thus, the P-wave energy gain 
is parametrically smaller in weak coupling. 
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Discussion. Several asymmetric phases proposed in 
the hterature contain Fermi surfaces, including the nor- 
mal Fermi liquid phases as well as the gapless breached 
pair phases. Kohn and Luttinger pointed out that, at 
zero temperature, all Fermi surfaces are unstable in the 
presence of interactions. We have shown that, in weak 
coupling, induced interactions lead to the formation of 
P-wave superfluids with maximal gaps for intermediate 
asymmetries. Thus, the suggested normal Fermi-liquid 
phases and breached pair phases are replaced by P-wave 
superfluids at zero temperature. In the asymptotic re- 
gions (both in asymmetry and coupling), the gaps are 
exponentially suppressed. The thermodynamic proper- 
ties of these phases will thus still be dominated by their 
underlying Fermi-liquid/breached pair nature. However, 
the dependence of the P-wave gaps on the S-wave scatter- 
ing length suggests that the induced P-wave superfluids 
may be observable as one approaches the Feshbach reso- 
nance, in particular from the BEC regime. 

The nature of both phases Pi_2 is that of P-wave pair- 
ing between fermions of the same species with order pa- 
rameter (a^a„g) Ap(p) = I]m^™^i:™(^p)- There 
are two qualitatively different states: one with two nodal 
points (Ap ^ ^L.ii) that breaks time-reversal invariance, 
and one with a nodal plane (Ap ^ ^1.0)- Energetic calcu- 
lations s ugg est that the state with broken time-reversal is 
favored [22|. Both states spontaneously break rotational 
invariance, so the low-energy spectrum will include a new 
Nambu-Goldstone boson |2J| in addition to the gapless 
nodes. In the P2 phase, there may be interesting en- 
trainment physics, since the two P-wave superfluids are 
coupled by the S-wave interaction. 

As shown in Fig.^ the phases Pi^2 will occupy shells 
in an optical trap at sufficiently low temperatures. As 
long as there is a finite molecular component rib, a thin 
shell of the Pi phase will always be realized. (The phys- 
ical volume may be maximized by tuning the densities.) 
Observations of Pi may also be facilitated by its occur- 
rence in the BEC regime. The P2 phase will occupy a 
large physical volume in the outer regions of the cloud 
containing an asymmetry in the BCS regime. One mech- 
anism for detecting these phases would be to look for 
rotational asymmetries in the momentum distribution or 
pairwise correlations |25l | in the time-of-flight expansion 
images of a single species. An interesting question is how 
the trap asymmetries would influence the spontaneously 
chosen orientation of the P-wave state. 
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